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We have calculated long-range molecular potentials of the Ojj~, 0^7 and l u symmetries between 
highly-excited rubidium atoms. Strong np + np potentials characterized by these symmetries are 
important in describing interaction-induced phenomena in the excitation spectra of high np Rydberg 
states. Long-range molecular resonances are such phenomena and they were first reported in S.M. 
Farooqi et al, Phys. Rev. Lett. 91 183002. One class of these resonances occurs at energies 
corresponding to excited atom pairs (n— l)d+ ns. Such resonances are attributed to ^-mixing due 
to Rydberg-Rydberg interactions so that otherwise forbidden molecular transitions become allowed. 
We calculate molecular potentials in Hund's case (c), use them to find the resonance lineshape and 
compare to experimental results. 
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I. INTRODUCTION 



Rydberg atoms have long been studied for their unique 
properties, such as long radiative lifetimes, large cross 
sections and huge polarizabilities At high principal 
quantum numbers, interaction forces between Rydberg 
atoms become extremely large. Manifestations of these 
long-range interactions can be seen in ultracold collisions 
or density-dependent line broadening of resonances in 
atomic beams j3fl . In ultracold Rydberg systems thermal 
motion is greatly reduced so that the effects of strong 
interactions can be investigated and utilized in a more 
controlablc manner. The prospect of using ultracold 
Rydberg atoms in quantum information applications has 
stimulated a great deal of experimental and theoretical 
interest lately [1, [f| . It was proposed to use the effect of 
dipole blockade to realize scalable quantum gates |5j. In 
mesoscopic ensembles of atoms, one excited atom pre- 
vents the Rydberg excitation of its neighbors because 
strong interactions shift many-atom excited states out 
of resonance. Recently, large inhibition of Rydberg ex- 
citation due to van der Waals (vdW) interactions (for 
n ~ 70-80), have been observed using apulse-amplificd 
single- mode laser @ and cw excitation [3, Q. Also, the 
dipole-blockade of cw Rydberg excitation of Cs atoms 
has been achieved [9(. 

There have been several proposals for weakly bound 
states involving Rydberg atoms but they have not yet 
been detected. The so called "trilobite" and "butterfly" 
states are molecular states formed by a pair of atoms, 
with one of the atoms in the ground state and the other 
one in a Rydberg state [13, El, E 0, El- Bound 
states of two Rydberg atoms have been also proposed 
[l5| . The long-range molecular resonances are another 
effect of Rydberg interactions. In the experiment [l6l |. 
(ri — l)p + (n + l)p and (n — l)d + ns resonances have 
been observed in single photon UV excitation from the 
5s ground state to high np Rydberg states (n = 50 — 90). 
A detailed theoretical treatment of (n — l)p + (n + l)p 
resonances was presented in pjj ■ However, the energy 



separations betwen the [n— l)d + ns resonances and the 
atomic np resonances are several times greater than the 
corresponding energy spacings for the (n— l)p+(n+l)p res- 
onances. As a result, much higher laser power is needed 
to excite the (n— l)d + ns resonances. Also, many more 
molecular states are needed for a detailed analysis of the 
interaction-induced ^-mixing, which is necessary for the 
existence of the resonances. 



II. THEORY 

This paper consists mainly of two parts. In the first 
part we calculate long-range molecular potentials of the 
0+, 0~ and 1„ symmetries since these symetries de- 
scribe strong np + np potentials. In the second part 
we evaluate the excitation dynamics of these molecular 
states to investigate the variouseffects of interactions be- 
tween the atoms. In this paper we primarily study the 
(n — l)d + ns resonances reported in [l(|. They occur at 
the average energy of excited atom pairs (n — l)d and 
ns, and do not correspond to any single- atom transi- 
tions. The contributions of diatomic potentials which 
coincide with the asymptotic (n — l)d + ns levels are 
dominant, although, there are significant contributions 
from (n — 2)d + (n + l)s and (n — 2)p\/2 + [n + 2)p 1 / 2 
potentials, with approximately the same asymptotic en- 
ergies. In one-photon transitions from the 5s ground 
state, dipole transitions to nd and ns states are not al- 
lowed. However, at high principal quantum numbers, 
long-range Rydberg-Rydberg interactions cause ^-mixing 
so that other molecular states, besides np + np, become 
accessible. Although the physics of this resonance and 
the (n — l)p+ (n + l)p one is very similar, the treatment 
of {n — l)d + ns resonances is technically much more 
demanding. There are many more asymptotic states be- 
tween np + np and (n — l)d + ns asymptotes and they 
all have to be included in order to describe the ^-mixing 
correctly. To make the potentials accurate at short dis- 
tances, many nearby asymptotes must also be included 
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in the asymptotic basis. In addition, the laser intensity 
used to excite them was almost two order of magnitude 
greater than that used for (n — l)p + (n + l)p resonances 
so there could be more power-dependent terms to con- 
sider, besides the two-photon Rabi frequency. We have 
used several different approaches to evaluate the results 
of the Rydberg excitation of atom pairs and they all give 
consistent results. 



A. Long-Range Rydberg-Rydberg potential curves 

We calculate long-range potentials in Hund's case (c) 
by diagonalization of an interaction matrix. There are 
several reasons why this Hund's case is considered. Here 
we point out only the basic arguments since they are 
explained in detail in [13]. We consider the states that 
are directly or indirectly coupled to np + np asymptotes, 
which can be directly excited from the 5s + 5s ground 
state. We can focus on strong potentials only because 
they can significantly mix, at short internuclear separa- 
tions, with other asymptotic states. Such np + np poten- 
tials are those coupled to nearby ns + (n + l)s states. 
The asymptotic energy spacing between np + np and 
ns + (n+ l)s states is quite small, which additionally in- 
creases the effect of their strong dipole-dipole coupling. If 
there were no atomic fine structure, the only candidates 
for strong np + np potentials would be the asymptotically 
degenerate 1 S+ and 3 £+ states [7]. After adding spin, 
the possible symmetries with strong potentials are 0^ , 0~ 
and l u . Fine structure can be, at least initially, neglected 
only if its energy splitting is much less than the separa- 
tion between adjacent nt + n'£' asymptotic levels, which 
is not the case here. For instance, the asymptotic energy 
spacing between 70p 1 / 2 +70p 1 / 2 and 70p 3 / 2 +70p 3 / 2 is al- 
most three times larger then the asymptotic separation 
between 70si/ 2 +71si/ 2 and 70p 3 / 2 +70p 3 / 2 . This obstacle 
impedes a unique definition of the dispersion coefficient 
Cq. However, these dispersion coefficients cannot be used 
because the interaction energy of pairs of Rydberg atoms 
contributing to (n— l)d + ns resonances is about twenty 
times greater than the region of energy for which pertur- 
bation theory is applicable. To obtain potential curves 
one has to diagonalize simultaneously both interaction 
terms, the long-range Rydberg term and the atomic fine- 
structure one. 

The interaction matrix includes long-range Rydberg- 
Rydberg interactions VR y d(i?) and the atomic spin-orbit 
interaction Hf s 



U(R) = V Ryd (R) + H i 



Is ■ 



(1) 



The eigenproblem of this interaction matrix is solved 
for the 0+, 0~ and l u states only. The first step is 
to construct asymptotic states of these symmetries and 
used them as a basis to represent U(R). The projection 
il = mi + rn 2 of the total angular momentum onto the 
molecular axis is conserved and properly symmetrized 



asymptotic states are constructed as follows 
\ntj,m,j;n'tj,,Q - rrij;Cl g / u ) 
~ \n,l,j,mj)\n\e\j'^-mj) 



(2) 



No overlap of the charge distributions of different atoms 
is assumed here since ultracold gases are very dilute, even 
considering these Rydberg atoms. All molecular states 
are defined in the molecule-fixed reference frame. States 
with Q = have to be (anti)symmetric under the action 
of the reflection operator <j v . We have addopted the fol- 
lowing convention [l8| for the action of a v 



<7„|A) = (-1) A |-A) , 
a u \S,M s ) = (-l) s - M s\S, 



-Ms 



(3) 
(4) 



The first rule ([3]) is obviously not applicable if A = 0. 
The correct result of <j„ on such states follows from its 
action on atomic states a v \£,m) = (— l) m \£, —m). 

After adding spin, there are many more different sym- 
metries to consider. On the other hand, it makes the 
^-mixing a smoother function of interaction energy, sim- 
ply because there are many more potential curves with 
a given ^-component for a given range of energy. We 
are primarily interested in how £ — 1 or p-character is 
distributed over different potential curves for each R. If 
a(R; A) is the fraction of a given asymptotic state in 
the molecular state |<p_r ; a), then the sum of a(R; A) over 
all molecular states has to be unity for each R. This 
statement is just the normalization condition of \<p) in 
the eigenbasis of U(R). We use this a(R;X) as a mea- 
sure of the character associated with \<fi). Therefore, for 
any np + np asymptotic state, the sum of the p-character 
of all molecular states has to be unity for any R. 

We have assumed in this analysis that there is no back- 
ground electric field. If a background electric field is in- 
cluded, then the group symmetry of the electronic Hamil- 
tonian is not -Dooh- Consequently, the quantum numbers 
based on this symmetry are not good anymore and we 
cannot reduce the interaction matrix in the basis of states 
of a given symmetry. In general, one has to include all 
possible basis states and use them to diagonalize the in- 
teraction matrix [l9| . 

For non-overlapping charge distributions, the 
Rydberg-Rydberg interaction may be expanded in 
a series of inverse powers of nuclear separations R. 
In this paper, we consider dipolar and quadrupolar 
interactions only. These interactions can be described 
by the following form [H U 



V L (R) = - 



(-r/47rrf rf 
LR 2L+1 



(5) 



where L = 1 (L = 2) for dipolar (quadrupolar) interac- 
tions, = m!/(n — m)! is the binomial coefficient, fj 
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FIG. 1: (a) Potentials curves for the 0^" symmetry. We 
present all the potential curves between 70s + 71s and 68s + 
73s. Molecular states, besides those that coincide with the 
asymptotic 70p + 70p states, may become accessible due to £- 
mixing induced by interactions. The resonance is dominantly 
produced by the 69d + 70s states, although significant con- 
tributions also come from the 68d -1- 71s and 68pi/2 + 72pi/ 2 
states. There are two 69d + 70s and 68d + 71s states of this 
symmetry and only one from the 68pi/2 + 72pi/ 2 asymptote. 
The fine structure of 69d and 68d states is ignored, (b) Av- 
erage radial dependencies of the two-photon Rabi frequency. 
The excitation probability of an atom pair depends on its sep- 
aration R since the fraction of p-character of any molecular 
state is a function of R. The radial dependence is shown for all 
states of the three asymptotes contributing to the resonance. 



is the position of the electron i around its center and 
L = 2L + 1. 

Defining \a)\b) = \n a ,£ a ,j a ,m a )\n b ,£ b ,j b ,il-m a ), one 
can show that the matrix elements of Vp,{R) are 



L-i-n+Y. n I 

<i,2|Vb|3,4> = (— l) «=* \jtxWmhhH 

„ ^n^ii (h L £ 3 \ f£ 2 L £ 4 \ \ j x L j 3 



R 2L+i [o o o A o o o J \ 4 \ £i 



32 L J 4 

£4 \ £2 



L+m 



h L j 3 

-mi m m 3 



32 



.74 



-fi + mi — m 57 — 7713 



{i\r L \j) 



where li = 2^ + 1, j t = 2j l 
part of the matrix element 

The asymptotic basis used to represent U (R) consists 
of all np+np, (n—l)d+ns and (n— 2)d+(n+l)s asymp- 
totic states, as well as all the states in between with 
a significant coupling to these asymptotes. These in- 
clude the (n — l)p + (n + l)p, (n — l)s + (n + 2)s and 
(n— 2)p+(n + 2)p asymptotes. These states are our pri- 
mary interest, but to describe them correctly at short in- 
ternuclear distances one should also include other nearby 
states strongly coupled to them. Included nearby asymp- 
totes are ns+(n+ l)s, (n — 2)s+(n + 3)s, (n—3)d+(n+2)s, 
nd+(n—l)s, (n— 3)f+np and (n— 2)/+(n— l)p. While the 
dipole-dipole interaction couples states belonging to dif- 
ferent n£ + n'£' asymptotes, the quadrupole interaction 
is mainly relevant for states within the same n£ + n'£' 
asymptote. The only exceptions are the off-diagonal 
quadrupole matrix elements between (n— l)d + ns and 
(n — 2)d + (n + l)s asymptotes. As mentioned before, 
the coupling depends on |ni — 7i2|, or to be precise, 
\ n i ~ n 2lj where the quantum defect 5? is included in 
the effective principal quantum number n* as follows: 
n* = rii — Se- Since the difference in effective principal 
quantum numbers of states (n— l)d and (n+l)s, as well 
as (n — 2)d and ns states, is only 0.22 for Rb (fo r high 
Rydberg states Si = q « 3.13 and 61=2 ~ 1-35 [22||), the 
off-diagonal quadrupole matrix element is several times 
larger than the diagonal ones. These asymptotes are very 
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, (6) 



Separation R (x 10 a.u.) 

FIG. 2: Same as Fig. [T]but for the 0^ symmetry 
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close in energy (separated by only 200 MHz for n = 70) 
so that at R ~ 30000 a and n = 70, quadrupole off- 
diagonal coupling is comparable with the asymptotic en- 
ergy spacing and these states become well mixed. This 
off-diagonal quadrupole coupling is relevant for the shape 
of (n— l)d+ns and (n— 2)d+(n+l)s resonances. 



B. Excitation of a pair of interacting Rydberg 
atoms 

We treat the (n— l)d+ns molecular resonances as two- 
body phenomena, which is supported by the fact that 
their energies coincide with the average value of only two 
atomic energies fl(| . The resonances are far red-detuned, 
so normally there should not be many excited atoms and 
presumably pair-wise excitation is the dominant mecha- 
nism. The general problem to be solved is the same as 
in [17j for (n — l)p + (n + l)p resonances. However, the 
analysis of (n — l)d + ns resonances is technically much 
more demanding and some approximations we have ap- 
plied before may not be satisfactory here. To check this, 
we have used several different ways to evaluate the con- 
tributions to the number of excited atoms from various 
molecular states. We discuss them in the results section. 

We consider a two-body Hamiltonian that includes 



70S+71S 
-70p+70p 




Separation R (x 10 a.u.) 

FIG. 3: Same as Fig. [T]but for the l u symmetry, except that 
each nd + n's asymptote has four states of this symmetry. 



long-range interactions and a linearly polarized optical 
field. For simplicity, we include only one molecular state 
| ip\) and the corresponding molecular potential e\{R). 
The Hamiltonian is (in the rotating frame and in the 
rotating-wave approximation) 



H 



2 

E 

i=l 



[A<7, 

+ [A A + eA (i?)]|^ A )(^ A 



2 

E 

1=1 



to 

r 



< y 



'e'g 



(7) 



where u>, lu' , and A, A' are single-photon Rabi fre- 
quencies and detunings relative to the nps/2 and npi/2 
fine-structure components, respectively. Here A A is the 
two-photon detuning from the asymptotic e\{R — > oo) 
molecular level. The operators a % and a\ e are de- 



a 'eg = J2 m \ e ^ m )(9i,m\ and a\ e 
To calculate lu and to' we use the ex- 

and 

2.3, which reflects 



fined as follows: 

Em \e l ,m)(e i ,m\ 

perimental values for the oscillator strength / 3 / 2 [2 
lo/lo' 



the ratio [24J w/u 1 = % // 3 / 2 //i/2 
the non-statistical character of the oscillator strengths 

h/2 and A/2- 

The excitation process is essentially a three-level 
scheme, although the number of states involved in the 
process is greater than three. We assume that a pair 
of ultracold Rb atoms is initially in one of the 5s + 5s 
ground states. There are four possible ground states 
1 5s, mj)|5s, "m'j), corresponding to different projections of 
spin nij — ±1/2. Ultimately, the total probability is av- 
eraged over all possible initial states. In the first excita- 
tion step, one of the atoms is excited to a given npj state. 
There are two intermediate states if the ground-state 
atoms have the same projections rrij, otherwise there are 
four of them. These intermediate states are further ex- 
cited in the second step. The final state in this excitation 
scheme can be a single molecular state \ip\(R)} or a su- 
perposition of states. We consider both cases but in the 
equations a single molecular state is assumed. Includ- 
ing more states in the equations is very straightforward 
and there is no need to do it explicitly. Eventually, to 
get the final excitation probability per atom for a given 
optical frequency, the contributions from all molecular 
potentials and all atom pairs that include a given atom 
are collected. 

To get pair excitation probabilities we solve the time- 
dependent Schrodinger equation for the ground diatomic 
state, all intermediate states and a given doubly-excited 
molecular state |</5a)- Here we present in detail the treat- 
ment if two rrij of the diatomic ground state are different 
and give only the final result when they are the same 
(this case is considered in detail in 17]). Utilizing the 
fact that symmetric and antisymmetric states have in- 
dependent time evolutions, the first step is to construct 
symmetric and antisymmetric combinations 



1 



\ij) = -j={\i,m) \ j,-m)+q\i,-m) \ j,m)}. 



(8) 



where q — 1(— 1) for symmetric (antisymmetric) states 
and m = 1/2 for linear laser polarization. In this way 
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the diatomic ground state \gg), four intermediate states 
\ge), \eg), \ge'}, \e'g), and four doubly-excited states |ee), 
|ee'), |e'e), |e'e') are denned. Here e and e' refer to ry>3/ 2 
and ropi/2 states, respectively. If there were no interac- 
tions, diatomic states |ee), |ee'), |e'e), |e'e') would be the 
states directly accessible in two-photon excitation. Any 
molecular state \ip\) is accessible if it has some compo- 
nents of these diatomic states. Due to ^-mixing induced 
by interactions, many \ip\(R)} gain a significant np frac- 
tion at some finite internuclear separation R. 
The wave function is modeled as follows 



=c \gg) + en \ge) + c 12 \eg) + c' n \ge') 
+ c 'i2 \e'g) + c 2 \(px) ■ 



(9) 



Solving the time-dependent Schrodinger equation 
idip/dt — H%\) (h = 1) leads to the coupled system for 
the excitation amplitudes c(t), 



(10) 



■ dc w* , lu'* . , 

den UJ uj* lu'* . 

«-^ = Acn + -c + Y< ee I^A)c 2 + — (e e\ipx)C2, (11) 

dc\2 . LU OJ* Lu'* .. . 

«-^ = Aci 2 + -c + y (ee|< / 3 A )c 2 + — (ee \ipx)c 2 , (12) 

J / / * /* 

i-^ = Ac , 11 + -c + — (ee> A )c 2 + — (e' e > A )c 2l (13) 
~dt 



*^ = Ac i2 + | c o + Y< e e|^A)c 2 + ^ r (eeV A )c 2 , (14) 



i^ = (Ax+e x (R))c 2 + (|<¥>A|ee) + ^ A | e ' e )j Cu 
+ (^(<P\\ee) + j (^|ee'))c 12 + (|<^ A |ee') (15) 
+ y(¥'A|e'e , ))c , 11 + (|<</? A |ee) + ^-(^ A | ee ')) c' 12 . 



The analogous system of equations for rrij = 



1/2 



was considered in detail [T3| in the analysis of (n— l)p 
(n + l)p resonances. The projections onto the molecular 
state are defined as: a ee (A) = (ee\tp\), a ee >(\) = (ee'\ip x ), 
fle'e(A) = (e'e\ip\), and a e > e >(\) — (e'e'\tp\). We assume 
that all of the ay coefficients are real. 

One can obtain a tractable system after the elimination 
of the excitation amplitudes of all intermediate states. 
This is justified by the fact that the dominant frequen- 
cies governing their time evolutions are A and A' which 
are much larger than the relevant Rabi frequencies. On 
a molecular resonance, A and A' are about 2tt ■ 2.2 GHz, 
while the peak values of u> and lu' are about 2.3 GHz 
and 1 GHz, respectively (for the actual experimental pa- 
rameters). Laser intensities used in this experiment were 
almost two order of magnitude higher than the ones used 
for the (n — l)p + (n + l)p resonances so that we have 
to check for the effects of higher laser fields. We adia- 
batically eliminate en, ci 2 , c' n and c' 12 , as in [llj, but 
this time keeping power-dependent terms. The result is 



equivalent to the Bloch equations of a two-level system 
.dc 



dt 
.dc 2 



-UJ g (t)CQ 



-c 2 , 



[A x +e x (R)-Lu e (t)]c 2 



Weff 



"Co, 



where the effective two-photon Rabi frequency is 

cu 2 lulu'A + A' 

^ c ff =-^a ee {\) 



(16) 
(17) 



AA 



— (a ee >(\) + a e / e (A)), 



1 2 



(18) 



•(A) 



and the power-dependent terms lu 9 and uj e are 



,'|2 



9 2A 2A' ' 

_ \ma ee + uj'a e ' e \ 2 
Je ~ 4A 

\to'a e i e > + Lua e 



4A 

\u'a e > e 



4A' 



4A' 



(19) 



(20) 



These tu g (t) and tu e {t) can be interpreted as power depen- 
dent shifts of the diatomic ground and doubly-excited 
states, respectively. Also, lu q is much greater than u) e 
because it does not depend on , which measure the p- 
character of the doubly-excite state. For the experimen- 
tal parameters, the peak value of ui g {t) on a molecular 
resonance is about 280 MHz, which is comparable with 
the laser bandwidth. 

In the vicinity of the molecular resonance, the second 
term on the right-hand side of Eq. (TT5]) is much smaller 
than the first one, and thus can be ignored. After ne- 
glecting that term, Eq. (fTB|) can be solved 



c (i) = exp 



i / u g (t')dt' 



(21) 



This Co can be used to find c 2 after the phase trans- 
formation c 2 = exp[— i(A\ + e\(R))t — J t uj e (t r ) dt']c 2 is 
performed. The excitation amplitude c 2 is 



c 2 



i l dti 

1 1 



w G ff(£i 



(22) 



x e i(A A + eA (iQ)ti+/4 1 (w g (i 2 )-Mi2))di2_ 



We assume that excitation laser pulses have a Gaus- 
sian time profile. It has been shown that the excess 
bandwidth in this experimental setup is mainly due to 
a linear frequency chirp. Because of the constant ratio 
lu/lu', both lu and lu' have the same time dependence 
lu,lu' ~ exp(— (1 + i~f)t 2 /a 2 , where a measures the pulse 
duration and 7 is a chirp parameter related to the laser 
bandwith T and duration a as follows 



ir 2 T 2 a 2 - 2 In 2 



2 In 2 



(23) 



6 



The probability to excite \<p\(R)) from the initial state 
\5s,m 1 / 2 }\5s,m 1 / 2 ) is -Pi (A) = \c 2 (t — > oo)| 2 . According 



to Eqs. 



and uj e do not depend on the chirp. 



If Lo g and uj e in Eq. (|22[) can be ignored, a simple formula 
can be derived [17| . 



2 _ /3 2 (A)tt J 2 ^ln 2 2 o „ [AA+£A(i?)1 2 /2;r 2 r 2 
1 21 2 J 2 r3r 



(24) 



where 7 sa t is the saturation intensity for ideal unchirped 
light and isolated np^/ 2 atoms and r is a Gaussian pulse 
duration (FWHM). The approximation applied to get 
this formula is equivalent to neglecting the w e -term in 
Eq. (jTTJ) and both terms in the right-hand side of Eq. 
(fl"tj)) (i.e. using c = 1). In Eq. ([24j), /3(A) is a time- 
independent part of w e s defined via Lu c g(t) ~ (3{X)uj 2 {t). 

If all a-cocfficients are equal to zero, according to Eq. 
(UHJ), the effective two-photon Rabi frequency is also 
zero. These coefficients measure different p-characters 
in \if\(R)). To evaluate them, we have to express the 
doubly-excited states |ee'), |ee), and |e'e') (defined in the 
space-fixed frame) in the molecule-fixed reference frame, 
where all \lp\(R)} are naturally defined [T3|- Apparently, 
all angular dependence related to different orientations 
of the molecular axis is contained in these a-cocfficients, 
and the observable quantities have to be averaged over 
all spatial orientations of the internuclcar axis. 

For the case nij = m'j = ±1/2, the ground and 
doubly-excited states are respectively \gg) = \g,m)\g,m) 
and |ee) = \e,m)\e,m) . Intermediate states \ge), \ge') 
and doubly-excited state |ee') have the form \ij) — 
-4= m)\j, m) + \j, m)\i, to)}. The coefficients a,j are 
defined as before. After the elimination of excitation 
amplitudes of all intermediate states, one gets again the 
system (|16H17|) with different uj c s(t) and to e (t) 



Wcff 



W , U)" ... 
— a ee (A)+— a e 'e' (A) + 



V2 AA 



-<w(A), (25) 



\\/2oja ee + uj'a ee i | 2 \y/2u'a e >e> + wa ee ' | z 



4A 



4A' 



(26) 



The solutions for these initial states are also given by 
Eq. (|22|) . or for low laser intensity by Eq. (f2"4|) . These 
excitation probabilities must also be averaged over all 
spatial orientations of the internuclear axis. In (|24| this 
is applied to /3(A) only. 

For an unpolarized sample of ultracold atoms, all ini- 
tial states are equally probable, so that excitation proba- 
bilities have to be averaged over all initial diatomic states 
and then all contributions to a given \tp\(R)} are summed 
up. This is done for many \<p\) to include the possibility 
of exciting different molecular potentials. The result is 
the averaged excitation P 2 (R) of an atom pair. The ex- 
citation probability per atom is the sum of all excitation 
probabilities of atom pairs that include a given atom, 



P fi xr. = 4tt 



dRR 2 pP 2 (R), 



(27) 



where p is the sample density. 

Molecular potentials can be very close to each other for 
some R, as shown in Figs. 1-3. For such R atom pairs 
are excited into superpositions of molecular states. Al- 
tough the extension of the system (|10j|15|) to include more 
molecular states \cp\) is simple, the technical difficulties 



^0.6 




2ph mol 




-2.6 



-2.4 -2.2 -2 

Detuning from 70p 3/2 (GHz) 



FIG. 4: (a) Fourier transforms of u>2 P h(i) of a chirped Gaus- 
sian pulse (dashed line) and its counterpart square pulse (solid 
line). The pulse areas and T2 P h (FWHM) of u>2 P h(i) for both 
pulses are matched to minimize the difference in the excita- 
tion probabilities. The tail of the square pulse vanishes much 
slower so it is truncated at A2 P h/r2 P h = 1.6, where A2 P h 
is the two-photon detuning. The remaining small difference 
between the pulses is expected to diminish further after sum- 
ming over all atom pairs in (|35[1 . (b) Comparison between 
calculated molecular signals using Eq. (|24[) (dashed line) and 
Eq. (|30|l (solid line) . The latter is based on the exact solution 
for a square pulse assuming that, in general, superpositions 
of molecular states are excited. The differences are somewhat 
larger just above the molecular resonance, where there are 
many very close molecular potentials, as shown in Figs. IH-'il 
The assumed laser bandwidth in this plot is 200 MHz. 
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of solving it are not negligible. Not only do we have four 
parameters to vary (i?, J, laser frequency and the ori- 
entation of the molecular axis), but we also have many 
molecular potentials and several symmetry cases. Also, 
for a given optical frequency, different molecular states 
are excited at different distances R. This is especially 
true at short R for which the potentials vary significantly. 
Instead of doing this, we choose a different approach to 
account for possible superpositions of molecular states. 
This approach is in many ways as complete as the full 
numerical calculation that includes the full set of doubly- 
excited states and at the same time is no more difficult 
than the calculation of molecular potentials itself. As Eq. 
(|24|) suggests, only a few parameters related to the exci- 
tation laser, such as the bandwidth and pulse duration, 
are important. The details of the laser pulse cannot be 
of fundamental importance, so we can substitute for the 
actual chirped Gaussian pulse a square pulse with the pa- 
rameters chosen to give probabilities consistent with (fM| . 
We have to match the pulse area and the width (actually 
FWHM) of the Fourier spectrum corresponding to the 
two-photon Rabi frequency w e ff. Therefore, V and r are 
not our direct concern, but rather r 2p h and T 2p h, which 
characterize Lo c s(t). For a Gaussian pulse T2 P h = V^r. 
This r 2p h would also be the FWHM of a square pulse if 
its duration were chosen to be T sq = 2.783/7rr 2 ph. The 
single-photon Rabi frequency uj sq of the square pulse is 
chosen to provide equal pulse areas of oj e g(t) of the ac- 
tual pulse and its substitute. A great adventage of this 
approach is that now we can easily write and calculate 
the exact excitation probabilities of all molecular states. 
The total Hamiltonian H to t of the system consists of the 
long-range interaction part U(R), given by Eq. ([1]) and 
the optical field (h = 1). 



H tot (R) =U(R) + ]T [A<4 + A'4, e ,] 

2 



s 1 ^ 

~2~ 69 



sq i 



e'g 



(28) 



To represent H to t we use the basis of U(R) completed by 
the intermediate states and the ground diatomic states. 
The matrix elements of H tot , corresponding to these 
added basis states, are essentially given by Eqs. (fTUl 
fT5|) . The only modification is to replace \ip\) by a su- 
perposition of different \<p\). The matrix of H to t is only 
slightly bigger than the matrix of U(R), so solving the 
eigenproblem of H to t does not impose additional difficul- 
ties. If £i(R) are the eigenvalues related to the eigenvec- 
tors \(j)i(R)) of H tot (R), then the solution of the time- 
dependent Schrodinger equation id \ ip) /dt = H tot \ip) is 



m;R))=Y,(MR)\99)e-^ (R)t \MR)) 



(29) 



Unlike Eqs. (|22l24p , the last formula does not give un- 
physical probabilities for large laser power. In the actual 



calculation we use a modification of Eq. ([29]) to account 
for the difference in excitation probabilities for large de- 
tunings. Excitation probabilities for a square pulse do 
not vanish sufficiently fast. This happens because the 
tails of the Fourier spectra of these two types of pulses 
are very different at large detunings. This is illustrated 
in Fig. Ufa). We overcome this by truncating the contri- 
butions from large detunings 



\ij e (t; r)) = Y^®ih(R) ' 2A| - r)T mol ) 

i 

x (MR)\gg)e- l ^ R)t \MR)) 



(30) 



where 8 is the Heavisidc function and rj is the cut-off 
parameter (our choice r\ — 1.6 is justified by Fig. Ufa)). 
We note that the radial dependence, shown in Fig. is 
obtained using the last equation. One can find the same 
dependence using the simpler formula (|22p . Both ways 
give essentially the same radial dependence. In the last 
formula all Si (R) are expressed with respect to the energy 
of asymptotic np 3 / 2 + np 3 / 2 states. 



C. n-scaling 

Now we estimate the major contribution to the n- 
scaling of the (n— l)d + ns resonance signal. For this 
purpose we use a rather simplified description of these 
states and the resonances. First we ignore fine struc- 
ture. As Fig. [5] suggests, the major contribution to the 
molecular signal does not come from the region of strong 
i'-mixing, but rather from the long-range region with con- 
siderably weaker mixing. In this estimate we completely 
neglect the contribution of strong ^-mixing at short R to 
the molecular signal. The wave function of (n— l)d + ns 
states can be expanded as follows 



\(n-l)dns;R) =\(n-l)dns {{ '>) + a pp (R)\npnp w ) 
+ \ip res ;R) , 



(31) 



where \(p res ; R) is the residue of the expansion 



and \ndn's 



'«(0)\ 



\npnp 



(0)\ 



are asymptotic states | (n 



l)dns; R— >oo) and \npnp; R^oo). In general, there are 
more than one \npnp^) state in the last expansion but 
the scaling law, in this approximation, does not depend 
of their number so keeping only one term is sufficient 
for our purpose. We want to find the function a pp (R) 
because the two-photon Rabi frequency is directly pro- 
portional to it. The Hamiltonian is still given by Eq. fl}. 
In the first approximation a pp (R) is 



o. 



(npnpW\V Ryd (R)\ndn'sW) 
Eds - Eppo 



(32) 



where Edso — Eppo is the asymptotic energy spacing of 
diatomic (n— l)d+ns and np+np levels. The last formula 
is valid in the region of weak £-mixing for (n— l)d + ns 
states. We conclude that 



,(R) ~ n 7 /R\ 



(33) 



8 



0.3 - 


1 1 


1 1 1 1 

/ \ Total 




i 1 

- 


' >- q 25 - 

o 
CD 

in 

o 0.2 ■ 


1 


La \ 




■ 


■i- 




\ \ 
' x \ 




■ 


| 0.15- 




i \ \ 
\ \ 
\ \ 




— 


Of O.I- 
'D ' 




\ \ 

1 \ \ 

/ " V 






0.05 ■ 


m 








n 1 


. As. .-O/T 


A 


~.~r.^.r_-_- 





2 4 6 8 10 12 14 



Separation R (x 10 5 a. u.) 

FIG. 5: Pair excitation probability as a function of sepa- 
ration R. We show contributions from all states associated 
with the three most relevant asymptotes. The total depen- 
dence includes twice the contribution of l u since this state 
is two- fold degenerate. According to Eq. (|27[l . dP cxc /dR = 
4/7r/9.R 2 P2(i?), where P2{R) is the excitation probability of an 
atom pair. Note that the weighted factor R 2 is included in the 
shown dependence. This figure shows that pairs at shorter 
distances are difficult to excite, even though I- mixing and 
their (u>1«) is much greater there (Figs. IH3[1 because they 
are very detuned from the resonance. The main contribution 
comes from pairs at separations R > 40000 ao- 

since only the dipole-dipole part of Vn y d(R) couples those 
asymptotic states, and Edso — Eppo ~ n -3 . 

We proceed using the results of the previous sections, 
but ignoring many details which are not of great impor- 
tance for n-scaling. The n-scaling is well denned only if 
there is no saturation of excitation so that the two-photon 
absorption probability Pi per pair is 

P2 ~ M 2 • (34) 

We have w e ff ~ <^ 2 app(-R)/A, where w is the single- 
atom Rabi frequency defined in the previous sections and 
A « (EdsQ — Eppo)/2 is the detuning from the atomic np 
resonance. Note that lo 2 /A is in the first approximation 
n-independent because the single photon Rabi frequency 
uo scales as the dipole matrix element n -3 / 2 and A ~ ri~ 3 . 
To get the excitation probability P cxc we use Eq. 

LO A f°° 

Pcxc^-rW dRR 2 \a pp (R)\ 2 ~n 14 /#o- (35) 

A JRo 

The lower limit Ro is set as follows. We assume that the 
laser frequency corresponds to the two-photon resonance, 
whose position coincides, to a very good approximation, 
with the asymptotic energy of the (n — l)d + ns state. 
A pair of atoms will be out of two-photon resonance if 
its interaction energy Ud s {Ro) is greater than Tiph ~ I\ 
therefore \Ud s {Ro)\ ~ E Because the laser bandwidth 
is considerably narrow (close to the Fourier transform 
limit), Ud s {Ro) is not very large, and assuming weak l- 
mixing, we can use second order perturbation theory to 



find an estimate for it. Eventhough the (n— 2)p+(n + 2)p 
states are very close to (n — l)d + ns states, they are 
very weakly coupled to them so they basically have no 
influence on this estimate. Ignoring the fine structure of 
(n— \)d+ns states, \Ud s (Ro)\ should be equal to \Cq\/R^, 
where Cq is the C coefficient for the (n — l)d + ns state. 
The estimate for the lower limit Rq is Rq 3 ~ yT/jQsj. 
Since Cg ~ n , we finally find 

Pcxc ~ n 8 ' 5 . (36) 



III. RESULTS AND DISCUSSION 
A. Numerical evaluation of excitation probabilities 

We have used four different ways to evaluate excitation 
probabilities in order to verify that the approximations 
are applicable for all conditions under which they are 
used. These conditions are different for different asymp- 
totes. One way is to use the method 17;] (based on 
Eq. (|23])) to get the lineshape of (n — l)d + ns reso- 
nances. However, it is assumed in this method that the 
effective two-photon Rabi frequency is sufficiently small 
and the molecular potentials are well separated. There- 
fore, for any pair of atoms at a certain separation R, 
only one doubly-excited molecular state is involved in 
the excitation (not necessarily the same one for all R). 
However, the two-photon Rabi frequency is almost an or- 
der of magnitude higher in this case so it is not obvious 
that some power-dependent terms can always be ignored. 
Such terms are included in Eq. f)22 f) . We find that, for 
most potentials, the formulae (ffi?)) - (|24[) give overall very 
similar lineshapes, with the position of the resonance 
slightly shifted, but the shape and amplitude are well 
preserved. However, it turns out that the approximate 
Eqs. (fM)) and ([^)) . for the experimental parameters, can- 
not be used for all laser frequencies and all asymptotes. 
By varying the laser frequency, we actually vary the re- 
gion of internuclear separations R for which atom pairs 
are on resonance. The two-photon Rabi frequency io c g is 
independent and, if the ^"-mixing is too large, the approx- 
imations may not be valid. To check if such parameters 
significantly affect the final probabilities, we have numer- 
ically solved the system (fTU)) - ([T5j) . 

It is solved for many different internuclear separations, 
spatial orientations of the molecular axis, laser intensi- 
ties and laser frequencies. This calculation is repeated for 
the states \ip\) which have dominant contributions to the 
resonance lineshape. They coincide with the asymptotic 
69c? + 70s, 68d+ 71s and 68pi/2 + 72pi/ 2 states. In these 
numerical calculations, only one molecular state \tp\) is 
considered as a final doubly-excited state of an atom pair. 
As explained previously, at some R, one should consider 
superpositions of molecular states \ip\). This is taken 
into account in Eqs. P5|) and (|3T))) . Even though we do 
not consider the actual pulse shape in this case, the pa- 
rameters of the substituted pulse are chosen to minimize 



9 



any quantitative difference between the pulses. Since we 
use exact solutions of 1(25]). this approach should give a 
rather fair description of the resonance phenomenon. 

It turns out that all these different methods lead to the 
same physical results, although, for some particular pa- 
rameters, they may be significantly different. Most of the 
differences vanish after summing over all potential curves 
and all atom pairs. The remaining variations are mainly 
due to the slight difference in the resonance positions and 
lineshapes obtained using different methods. The real 
physical parameters, such as the linewidth, signal size, 
rt-scaling or even the resonance position, are practically 
unchanged, as shown in Fig. |U Here we present the re- 
sult based on Eqs. (|2"9"j) - (|37)|) . the most complete method 
we use, and the results from the simplest method based 
on Eq. (j24|) . This gives an estimate of the variations 
between these methods. Our calculation shows that the 
probabilities of the 0~ asymptotes are about one order 
of magnitude less then the contribution of the other two 
symmetries. Interestingly, for (n — l)p + (n + l)p res- 
onances, the contribution of the 0+ symmetry was in- 
significant [171 ] . 



B. Comparison with theoretical lineshapes 

In Figs. QJa)-[3ja), we show three sets of molecular 
potentials corresponding to the three symmetries con- 
sidered. At short distances these potentials have very 
complicated shapes due to multiple avoided crossings. 
At these avoided crossings the ^-mixing is the strongest. 
The np+np components of various nearby potentials are 
the most important for the molecular resonance. The 
relevant physical quantity which depends on the frac- 
tion of np + np states is the two-photon Rabi frequency 
w c ft, defined by Eq. ([18]) . Asymptotes (n— l)d + ns, 
(n — 2)d + (n + l)s and (n — 2)pi/ 2 + (n + 2)p 1 / 2 give the 
essential contributions to the resonance. In parts (b) of 
Figs. [l]|3]we present the magnitude of the radial depen- 
dence w e ff for the three asymptotes and all their states. 
For each of these states, the radial dependence is ob- 
tained after averaging w^g over different orientations of 
the molecular axis and initial states. Even though lo c s for 
atom pairs at short distances is larger due to stronger l- 
mixing, such pairs are difficult to excite because they are 
very detuned from the molecular resonance. The actual 
pair excitation probability as a function of R, on exact 
molecular resonance, is shown in Fig. The weighting 
factor R 2 is included in the presented dependence. Sur- 
prisingly, there is a significant contribution to the molec- 
ular signal from the pairs at larger distances. 

We illustrate different methods used to calculate the 
resonance lineshape in Fig. 0Jb). We present the sim- 
plest method, given by Eq. (fM| . and the method based 
on the exact solution for square pulses (|29ti30[) . Only the 
latter method allows superpositions of molecular states 
to be excited. This possibility is relevant if potentials 
are very close to each other. These two methods give 
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FIG. 6: (a) Scaling of the position of molecular (n — l)d + ns 
resonances with respect to the atomic np 3 / 2 resonance. The 
position follows the characteristic n~ 3 scaling, (b) Scaling of 
the molecular signal. Here we test the expected n 8,5 scaling 
(I36p . This plot shows a reasonable agreement with the n 8 ' 5 
scaling. 



overall very similar lineshapes, and the difference is only 
in the details. The relative difference between them is 
somewhat larger just above the resonance, where there 
are many potentials very close to each other. For the 
position of the molecular resonance, the method method 
based on Eqs. (|29ti30[) gives —2.18 GHz from the atomic 
70p 3 / 2 resonance, which is in agreement with the experi- 
mental position of 2.21(3) GHz. It appears that the posi- 
tion of the resonance obtained using the simplest method 
is closer to the experimental value. However, this is likely 
just a coincidence because, in that calculation A and A' 
were replaced by their average value, so that fine details 
of the atomic level positions were not included. 

We have tested the n-scaling of the molecular reso- 
nance, both numerically and experimentally. As men- 
tioned, this scaling law makes sense only if, for a given 
laser intensity, the two-photon transition is not saturated 
for all values of n for which the scaling law is used. The 
calculated ratio (using Eq. ([24]) ) of the signals for all 
n = 70 and n = 60 is 3.915, while (70*/60*) 8 5 = 3.920. 
Also the same ratio for n = 90 and n = 70 is 9.274, 
while (90*/70*) 8 - 5 = 9.115. The experimental depen- 
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dence for the n-scaling of the molecular signal is shown 
in Fig. [B^b). The agreement is fairly good and the de- 
viation could be explained by variations of experimen- 
tal parameters between these three experimental scans. 
The n-scaling of the resonance position is in an excel- 
lent agreement with the expected scaling law, as shown 
in Fig. Ufa). 

In Fig. [7Ja) we present calculated lineshapes for sev- 
eral laser bandwidths, as indicated in the graph. The 
linewidth of the resonance is significantly larger than the 
laser bandwidth and is primarily determined by the de- 
tails of the molecular potentials and ^"-mixing. The con- 
voluted lineshapes for a laser bandwidth of 200 MHz is 
shown in Fig. El^b). We get the best agreement with the 
experiment for this bandwidth. The actual laser band- 
width was probably smaller than this one. For the ac- 
tual experimental conditions, this theory cannot be used 
to fit the portion of the red tail of the spectrum closer 
to the atomic resonance because the excitation fractions 
are much larger. The presence of excited atoms could 
modify the pair excitation probability so that Eq. (12"T1) 
is not applicable. Also, the simple pair excitation model 
cannot be used to explain the excitation process at larger 
excitation fractions. At such fractions, I- mixing could be 
more efficient because close molecular potentials could 
additionally mix due to the interactions with nearby ex- 
cited atoms. This could explain why the experimental 
red tail, for the experimental conditions, is just a simple 
monotonic function even though there are many poten- 
tial curves and avoided-crossings between the np+np and 
(n — l)d + ns asymptotic levels. 

We compare experimental and theoretical signal sizes 
assuming a typical 5 x 10 10 cm~ 3 density. The laser in- 
tensity is typically in the range of 400-500 MW cm~ 2 . 
Experimental signals are about 300 ions per shot. The 
details of the experimental setup are given in [l6| . From 
the scaling factor introduced in order to compare the 
theoretical and experimental lineshapes, we find that the 
calculated number of atoms is about 6-7 times higher 
than the experimental value. This is probably accept- 
able considering the number of factors which influence 
this estimate. Besides uncertainties in the experimen- 
tal parameters, this scaling factor may also suggest that 
there was some excitation blockade [6j, Ijj . 

IV. CONCLUSION 

We have presented long-range doubly-excited molecu- 
lar potentials of the 0+, 0~ and 1„ symmetries. These 
potentials are important in describing the effects of in- 
teractions in single-photon excitation to high np Rydberg 
states. We have illustrated the ^-mixing induced by in- 
teractions over a broad range of internuclear distances. 
Several methods to evaluate the spectrum of the excita- 
tion of interacting Rydberg pairs of atoms were presented 
and they all gave similar results. The analysis showed 



that molecular resonances at the average atomic nd and 




1/2 of Offset from 69d+70s resonance 

FIG. 7: (a) Calculated molecular signal for different laser 
bandwidths of the excitation laser, as indicated in the plot. 
The linewidth of the resonance is dominantly determined by 
the details of long-range interactions and it is significantly 
larger than the laser bandwidth, (b) Comparison of the the- 
oretical signals shown in Fig. |4jb) and the experimental one. 
This is the best fit of the experimental data, although the 
actual bandwidth was likely smaller than the one we used for 
the comparison. 

ns energies are expected to occur. The calculated prop- 
erties of these resonances, such as position, linewidth, 
n-scaling and signal size are reasonably close to the ex- 
perimental observation, but a complete understanding of 
the spectral features of Rydberg excitation requires im- 
proved theoretical models. 
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